Spontaneous trimerization in a bilinear-biquadratic S = 1 zig-zag chain 



Philippe Corboz/ Andreas M. Lauchli,^ Keisukc Totsuka,^ and Hirokazu Tsunctsugu'* 

' Institut fiir Theoretische Physik, ETH Zurich, CH-8093 Zurich, Switzerland 
^Institut Romand de Recherche Numerique en Physique des Materiaux (IRRMA), CH-1015 Lausanne, Switzerland 
' Yukawa Institute for Theoretical Physics, Kyoto University, 
Kitashirakawa Oiwake-Cho, Kyoto 606-8502, Japan 
''^Institute for Solid State Physics, University of Tokyo, Kashiwa, Chiba 277-8581, Japan 

(Dated: February 1, 2008) 

Recent theoretical studies raised the possibility of a realization of spin nematic states in the 5 = 1 
triangular lattice compound NiGa2S4. We study the bilinear-biquadratic spin 1 chain in a zig-zag 
geometry by means of the density matrix renormalization group method and exact diagonalization. 
We present the phase diagram focusing on antiferromagnetic interactions. Adjacent to the known 
Haldane-double Haldane and the extended critical phase with dominant spin nematic correlations 
we find a trimerized phase with a nonvanishing energy gap. We discuss results for different order 
parameters, energy gaps, correlation functions, and the central charge, and make connection to field 
theoretical predictions for the phase diagram. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.40.Cx, 75.40.Mg 



Introduction — Quantum spin systems have provided a 
very wide playground for the quest of novel quantum or- 
ders, and the short catalog includes Haldane gap, dimer 
order, chiral order and others. Recently the discovery of 
spin liquid like behavior in the spin-1 triangular magnet 
NiGa2S4 [1] has stimulated increasing interests in another 
type of quantum order: spin nematic order. This is the 
long range order of quadrupolc moments of local spins, 
in contrast to dipolc moment order parameter in con- 
ventional magnetic long range orders. In a spin nematic 
ordered phase, although spins show no static moment, 
spin rotation symmetry is spontaneously broken due to 
anisotropic spin fluctuations. For S=l spin operators, 
anisotropic spin fluctuations still have uniaxial symme- 
try, and this symmetry axis is called director, in analogy 
to a liquid crystal. Ferro and antiferro spin nematic or- 
ders were independently proposed by different groups as 
an explanation for the unusual low-temperature proper- 
ties of NiGa2S4 [2, 3, 4]. In particular, the antiferro spin 
nematic order is possible to match the triangular lattice 
structure without any frustration, and the ground state 
is unique aside from degeneracy due to global spin rota- 
tion. Although it remains open if antiferro spin nematic 
order is realized in this material, it is interesting and 
also important to investigate further this state and ob- 
tain a better understanding. To this end, we shall study 
a one-dimensional analog with the same three sublatticc 
structure, a zig-zag chain. 

Model — The Hamiltonian of the bilinear-biquadratic 
(BLBQ) spin 1 zig-zag chain is given by 

H ^ Ji ^cos6'(Si-Si+i)-hsin6l(Si-Si+i)^ 

i 

+ J2 ^cos0(Si-Si+2) + sin0(Si-Si+2)^ (1) 

i 

Si's are spin 1 operators and 6 parametrizes the strength 
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FIG. 1: Zig-zag geometry of the model with Ji (J2) the near- 
est (next nearest) neighbor coupling strength. 



of bilinear and biquadratic coupling. For 6 = {6 = 
7r/2) the biquadratic terms (bilinear terms) vanish. For 
6 ~ tt/A the Hamiltonian exhibits SU(3) spin symme- 
try. Ji and J2 are the nearest and next nearest neighbor 
coupling strengths, respectively. The model is best visu- 
alized in a zig-zag geometry where the Ji bonds couple 
two chains and the J2 bonds are located along the chains 
(Fig. 1). We concentrate on antiferromagnetic interac- 
tions on all bonds with < 9 < tt/2 and Ji, J2 > 0. By 
setting 9 = the model reduces to the Ji — J2 spin 1 
chain which exhibits a first order transition between two 
distinct topological orders [5]. The Haldane phase for 
J2/J1 < at ~ 0.775, which also contains the well studied 
spin 1 chain for J2 = 0, is characterized by a finite string 
order parameter (SOP). The ground state is a valence 
bond solid (VBS) where each site consists of two spin 
1/2 forming a triplet and two spin 1/2 on neighboring 
sites couple to a singlet. The transition into the double 
Haldane (DH) phase corresponds to a decoupHng of the 
VBS string into two intertwined substrings, where each 
string exhibits string order [5], leading to a finite double 
string order parameter (DSOP). 

The line J2/J1 = corresponds to the BLBQ chain 
studied, e.g., in [6]. The gapped Haldane phase (includ- 
ing the Affleck-Kennedy-Lieb-Tasaki (AKLT) point for 
tan 6* = 1/3) ends at the SU(3) symmetric point 9 = 7r/4 
(exactly solvable Uimin-Lai-Sutherland model) where a 
second order phase transition of Berezinskii-Kosterlitz- 
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FIG. 2: (Color online) The phase diagram of the spin 1 BLBQ 
zig-zag chain. The area of the black dots in the trimerized 
phase scales with the magnitude of the trimer order param- 
eter. We verified that the central charge is 2 in the critical 
phase for several points (red diamonds). The green triangles 
mark the phase boundary of the trimerized phase obtained 
from a level spectroscopy analysis from the ED data. 



Thouless (BKT) type occurs [7]. Thus there is an expo- 
nentially slow opening of the gap on the Haldane side. 
For 7r/4 < 6 < tt/2 one finds an extended critical phase 
with soft modes at fc = 0, ±27r/3 and central charge 
c = 2. The dominant correlations in this phase away from 
the SU(3) point are quadrupolar (spin nematic) [6, 7]. 

Numerical simulations — The numerical results have 
been obtained by the density matrix renormalization 
group (DMRG) [8, 9] method and exact diagonaliza- 
tion (ED). For the DMRG calculations we used up to 
m = 2000 states with typically 6 sweeps and system 
sizes with open boundary conditions up to L = 300 sites. 
Quantities of interest are extrapolated in m and if needed 
also for L — > oo. Error bars are estimated from the con- 
vergence behavior in m. With ED we considered systems 
with periodic boundary conditions up to L = 21. 

The overall phase diagram — Figure 2 summarizes 
the phase diagram of the Hamiltonian (1) obtained by 
DMRG and ED simulations. We reproduced the results 
for the phases on the axes: Haldane, the DH, and the ex- 
tended critical phase with c = 2. All these phases extend 
into the plane. Interestingly they all touch the dominant 
phase in the center, a gapped trimerized phase, which 
will be discussed below. 

The Haldane - double Haldane transition [5] point ex- 
tends as a first order line in the parameter space which 
terminates upon touching the boundary of the trimerized 
phase. We confirmed the first order nature of the tran- 
sitions by calculating the SOP and DSOP along several 
cuts for fixed 9 and varying J2/ Ji- 

Trimerized phase — The most exciting feature of the 
phase diagram is the gapped, trimerized phase, where 



three neighboring spins couple predominantly to a sin- 
glet. The trimer ground state is threefold degenerate 
and breaks translational invariance. This phase - includ- 
ing the frustration process leading to it - is reminiscent 
of the dimerized phase of the J1-J2 spin 1/2 chain for 
J2/J1 > 0.2411 [10]. 

Initially a massive trimerized phase for the spin 1 
Heisenberg chain (J2/J1 = 0) for 7r/4 < 9 < tt/2 was 
put forward [11, 12], but later works [6, 7, 13, 14, 15, 16] 
showed that the region remains massless and has domi- 
nant nematic correlations. In our model the additional 
next nearest neighbor coupling J2 allows us to stabilize 
the trimerized state. 

In previous work [17, 18] parent Hamiltonians have 
been constructed using complicated four site interactions, 
which exhibit exact trimer ground states (in Ref. [19] yet 
a different kind of trimerized state was constructed) . We 
now have found a trimerized phase that is stable in a 
finite region of the parameter space of a much simpler 
and possibly realistic spin Hamiltonian. We expect our 
model (1) to be related to these parent Hamiltonians in 
Refs. [17] and [18] in the same spirit as the Heisenberg 
spin 1 chain is related to the AKLT model. 

One of the most direct indications for a trimerized 
phase is the period 3 in the local bond energies (insets of 
Fig. 3). This pattern is formed because the two bonds be- 
longing to a trimer have a lower energy than the bonds 
connecting two trimers. We determine this oscillation 
amplitude of the local bond energies in the middle of the 
chain for different system sizes L and extrapolate it to 
L ^ 00. In the Haldane (respectively DH) phase the 
amplitude vanishes exponentially with L. In the critical 
region the amplitude extrapolates to zero with a power 
law. But in the trimerized phase the extrapolation of the 
oscillation amplitude clearly yields a finite value, which 
we call the trimer order parameter (TOP) (see Fig. 3). 

The magnitude of the TOP in the trimerized phase is 
proportional to the area of the black dots in Fig. 2. On 
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FIG. 3: (Color online) Extrapolation of the bond energy oscil- 
lation amplitude in the middle of the chain leading to a finite 
trimer order parameter in the trimerized phase (9 — 0.287r, 
J2/Ji=l). Inset: The local bond energies form a pattern of 
period 3 (L = 48 in this example). 
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the SU(3) line DMRG predicts that the trimerized phase 
sets in at J2/J1 « 0.45 and ends at « 3.5. We used 
a level spectroscopy analysis of the ED data to comple- 
ment the results (green triangles in Fig. 2). This tech- 
nique has been successfully applied in the case of the 
spin 1/2 chain to locate the transition point into the 
dimerized phase [10]. ED yields a consistent result with 
DMRG for the lower boundary of the trimerized phase, 
but it suggests a substantially smaller value for the upper 
boundary. This mismatch stems partly from strong finite 
size effects on the location of the level crossings (which 
arc much stronger than in the spin 1/2 case). We ob- 
serve that the ED phase boundary shifts to slightly larger 
J2 values upon taking bigger system sizes into account, 
approaching somewhat the DMRG results. A different 
source of error is the rather slow convergence of DMRG 
in m in or close to a critical phase, especially for large 
J1I J\- The important result, however, which we get from 
both methods is the finite extent of the trimerized phase 
on the J2/J1 axis, in contrast to the spin 1/2 case where 
the dimerized phase extends to infinity. We comment on 
a possible explanation in the field theory section below. 

Figure 4 shows the finite spin gap in the trimerized 
phase, /^E = Eoil) - So(0) where £^0(5"^) is the ground 
state energy in the sector. On the SU(3) line the 
lowest excitation is eightfold degenerate (i.e. spin and 
quadrupolar excitations are symmetry related on this 
line). Away from the SU(3) line the lowest excitations 
have Stot = I ioi 9 < 7r/4 and Stot = 2 for 6* > 7r/4. 
According to the analysis of Ref. 18, the nature of excita- 
tions in the trimerized phases are gapped, deconfined do- 
main walls, very similar to the deconfined spinous of the 
spontaneously dimerized phase of the frustrated 5=1/2 
spin chain. 

The BKT transition line lies on the SU(3) line except 
for 0.45 < J2/J1 < 3.5 where it follows the right bound- 
ary of the trimerized phase. Field theory predicts that 
the transition line separating the Haldane (respectively 
DH) phase from the trimerized phase is of first order 
(see below). We confirmed this numerically from a non- 
vanishing spin gap and a non diverging spin-spin corrcla- 
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FIG. 5; (Color online) Energy gap between the ground state 
of singlet and triplet sector (upper plot) and spin correlation 
length (lower plot) for different system sizes for fixed J2/J1 = 
1.25. The finite energy gap at the transition and the weak 
size dependence of the correlation length indicate a first order 
phase transition. The correlation length was obtained by a fit 
to the exponential decay of the spin-spin correlations. 



tion length across the transition for J2/J1 = 1.25 (Fig. 5). 
For J2/ Ji = 0.75 we cannot completely rule out a vanish- 
ing gap as it is rather small at the transition. Another 
indication for the first order nature is the jump of the 
TOP at the phase boundary, in contrast to the exponen- 
tial suppression of the TOP towards the critical phase. 

Critical phase — We have determined the value of the 
central charge numerically from the entanglement en- 
tropy for different points in the critical phase (red di- 
amonds in Fig. 2), and all the values agree with c = 2 
(to within 5%), corresponding to a level-1 SU(3) Wess- 
Zumino-Witten (WZW) model [7] . We made use of the 
universal scaling behavior of the entanglement entropy 
for conformally invariant one-dimensional quantum sys- 
tems [20] of the form S{x, L) = | In [^ sin(^)] + const, 
where S{x) is the von Neumann entropy of the reduced 
density matrix p{x) of a subsystem starting at the open 
boundary of length x embedded in a system of length 
L, and c is the central charge. In the critical phase 
the spin and quadrupolar correlation functions are ex- 
pected to decay as a power law with distance x with an 
exponent 77 = 4/3 and multiplicative logarithmic correc- 
tions [7] ^ cos(27ra;/3)(lna;)'^/a;'', which suppress the spin 
{a = —2) and enhance the quadrupolar (cr = 2) correla- 
tions, such that the latter are dominant (Fig. 6), similar 
to the single chain case [6] . 

Field theoretical considerations — As is discussed in [7, 
21], the low-energy effective Hamiltonian for the model 
(1) around the point J2 = 0, 6* = 7r/4 is the level-1 SU(3) 
WZW model perturbed by marginally irrelevant current- 
current interactions: 

nwzw+vgi JdxJ2 JlJr+v92 J dx J^J^ , (2) 



FIG. 4: (Color online) Energy gaps of spin excitations along 
cuts in the parameter space with J2/J1 = 1 (left-hand side) 
and 9 = tt/4 (right-hand side). 
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where the first summation is over the S0(3)-subset of the 
SU(3)-currents Jf^p. The choice ffi = 32 = 3 recovers 
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FIG. 6: (Color online) Spin (stars) and quadrupolar (crosses) 
correlation functions (normalized) for 9 — 0A2n, J2 = 0.25. 
The two solid lines are fits to the functional form predicted 
by field theory. Logarithmic corrections to the power law lead 
to dominant quadrupolar correlations. 



SU(3) symmetry and the intcgrablc model (J2 = 0,6 = 
7r/4) corresponds to g < 0. 

At the leading order, the effect of J2(>0) (with 9 = tt/A 
fixed) may be taken into aecomit by increasing g{=gi=g2) 
in Eq. (2). At some critical value of J2, g changes its 
sign to positive and the system flows into a massive phase 
with three degenerate ground states [22]. The model (2) 
has another massive phase characterized by the asymp- 
totic trajectory gi = —g2, which has a unique (Haldane) 
ground state and is separated from the trimerized one 
by a first-order transition [22]. The Y-shaped boundaries 
(among trimerized, Haldane and critical) and other quali- 
tative features in the lower portion of Fig. 2 are consistent 
with what is predicted by the model (2). The appearance 
of the gapless phase in the large- J2 region (J2/J1 > 3.5), 
as revealed by DMRG and ED may be understood as 
follows. Let us start with two weakly coupled BLBQ 
chains. From the field-theoretical viewpoint, the ordi- 
nary two-leg ladder and the zig-zag chain should behave 
similarly in the decoupled chain limit; both systems will 
fiow toward a strong-coupling limit upon switching on 
the interchain coupling Ji . From the standard expansion 
with respect to the strongly coupled rungs, we know that 
the strong-coupling phase of the ordinary two-leg BLBQ 
ladder {9 > 7r/4) is critical (as in the S = 1/2 three- leg 
ladder). Therefore we may expect the same critical be- 
havior in the zig-zag ladder as well at least for sufficiently 
large J2/J1. 

Conclusions — We have presented a study of a spin 
1 generalization of the famous Ji — J2 5'=l/2 model, 
motivated by recent proposals for spin nematic ground 
states in a spin 1 triangular lattice [2, 3, 4]. The phase 
diagram generalizes the well-known spin fiuid-dimerized 
transition to a (nematic) spin fiuid to trimerized tran- 
sition in the level-1 SU(3) WZW universality class. We 
further explored the phase diagram in the vicinity of this 
transition, revealing a realization of the gapped sector of 
the Andrei-Destri model. Coming from the limit of two 
decoupled BLBQ chains the interchain interaction is rele- 



vant and drives a crossover to gapless single-chain behav- 
ior, in contrast to the 5=1/2 case, where the marginal 
interaction lets the system flow to a dimerized strong 
coupling state. 

As a perspective we believe that the SU(N) Ji — J2 
model contains a general mechanism where the critical 
state realized for the Ji model flows towards a 7V-merized 
gapped state when the ratio J2/J1 is beyond a certain 
critical value. Another interesting point is the possibil- 
ity of stabilizing trimerized phases in two dimensional 
lattices, possibly in a Ji — J2 BLBQ model on the tri- 
angular lattice. Furthermore in the light of the present 
study, where we revealed dominant quadrupolar correla- 
tions in the vicinity of the SU(3) line, it will be interesting 
to explore whether other spin nematic or spin-multipolar 
phases [23] remain to be uncovered close to SU(N) regions 
in higher spin antiferromagnets. 
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